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Abstract 

> . 

, We show that gravitinos produced by decays of a supersymmetry breaking scalar field 

(the pseudo-moduli field) can naturally explain the observed abundance of dark matter 
in a certain class of the gauge mediation models. We study the decay processes as well 
as cosmological constraints on this scenario in detail, particularly focusing on different 
behavior of the real and imaginary components of the pseudo-moduli field. Cosmologically 
viable scenario emerges when the gravitino and the pseudo- moduli masses are 0(10 — 



100) MeV and 0(100) GeV, respectively. 
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1 Introduction 



The presence of dark matter (DM) in the Universe was firmly established by numerous 
observations jl]. Nevertheless, it has remained as a big mystery in cosmology as well as 
particle physics what DM is made of. Since there is no candidate for DM in the standard 
model, we need to consider new physics. 

In a supersymmetric extension of the standard model, the lightest super symmetric particle 
(LSP) is a stable particle if the R parity is conserved. Depending on the mediation mechanism 
of the supersymmetry (SUSY) breaking effects, the lightest neutralino or the gravitino can 
thus be a good candidate for DM. The latter possibility is naturally realized in a framework 
of gauge mediation |2|, [|, which has a virtue of avoiding the SUSY flavor problem. 

The production mechanisms of gravitinos in the early Universe are broadly classified into 
thermal or non-thermal one. The thermal production is always present as long as the Universe 
becomes radiation dominated after inflation Q ^, |(| [?], ||, || [l(J. In this case the decay rate of 
an inflaton must be such that gravitinos, produced by particle scatterings in thermal plasma, 
account for the observed DM abundance. If the inflationary dynamics has nothing to do with 
the SUSY breaking mechanism (and therefore the gravitino mass), such a coincidence may 
call for some explanation. On the other hand, non-thermal gravitino production has been 
discussed (mostly as a problem of overproduction or a solution to it) in the context of the 
decay of the next to lightest SUSY particle (NLSP) 0, Q , the moduli @ 0, [Tj| f[jj 0, [jlj , 



the inflaton fig, |20|, |2l|, |22|, £3|, |2J, |2§ g§ and the SUSY breaking field (sometimes called 



as the Polonyi field or the pseudo-moduli field) |27], |28|, |29|, |3C|, |3l|, [32|, |3J]. In particular, 
it is interesting to see if a right amount of the gravitinos can be produced by the decay of 
the SUSY breaking field, since the structure of the SUSY breaking sector may be probed by 
cosmological arguments. 

In this paper, we investigate the gravitino DM scenario in a generic setup of gauge- 
mediated SUSY breaking models. In many SUSY breaking models, there is a light singlet 
scalar field, which obtains a mass from SUSY breaking. During the inflation era, this scalar 
field, the pseudo-moduli field, can have a large displacement from the true vacuum, and 
at a later time, it starts coherent oscillations about the minimum of the potential. Under 
reasonable assumptions, the oscillation energy dominates the energy density of the Universe, 
and the decay of the scalar field produces radiation as well as gravitinos which remain as DM 
today. We update the calculation of Ref. |3^] by taking into account the following points. 
We do not assume a particular relation among parameters in the SUSY breaking sector. We 
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treat the real and the imaginary parts of the scalar field separately as their decay properties 
are quite different. We find that, for the mass of the scalar field around 0(100) GeV and the 
gravitino mass of O(10 — 100) MeV, the decay of the imaginary part is the main source of the 
radiation and the gravitinos that account for the observed DM abundance. The region turns 
out to be similar to the one found in Ref. |3^] although the main decay mode is different. The 
consistent region overlaps with the prediction of a model in Ref. |33| where the \i problem is 
solved .Q 



2 Gauge mediation 



We first define the framework and identify parameters relevant for the discussion of cosmology. 
We use an effective description of gauge- mediation models given in terms of a SUSY breaking 
field S and the fields in the minimal SUSY standard model (MSSM). The SUSY breaking 
sector is described by a single chiral superfield S which consists of the Goldstino fermion and 
its scalar partner s with the Kahler and super-potentials: 

K = S^S — ^ S ^J 2 + (higher order), (1) 

W = m 2 S, (2) 

where A is a cut-off scale of the effective theory, and m denotes the size of the SUSY breaking. 
This form is obtained after integrating out massive fields in a wide class of SUSY breaking 
models. The equation of motion gives F$ = m 2 as long as there is no singularity in the 
Kahler potential. The second term in Eq. (Q) stabilizes the scalar potential at s = 0. 

The MSSM particles can couple to the SUSY breaking sector through messenger fields, / 
and /: 

W 3 -XSff, (3) 

where A is a dimensionless coupling constant. With this term, the potential is minimized 
at s = and // = m 2 /X where SUSY is unbroken and the gauge symmetry of the MSSM 
is broken. Therefore, one needs some mechanism to stabilize the potential at (s) ^ and 
</) = (/) = o. 



*A similar parameter region is identified in Ref. J34j in the F-theory GUT model by considering the 
abundance of gravitino dark matter and the ^t-problem. In discussions of cosmology, the most important 
difference between two models is the mass of the imaginary part of the SUSY breaking scalar field. In the 
F-theory GUT model, the imaginary part is assumed to be much lighter than the real part. In the model of 
Refs. [B5L p2l B3|, in contrast, the real and the imaginary parts have almost the same masses. 



3 



In order to keep discussion as general as possible, we do not specify such a mechanism in 
the following and treat the three quantities ((s), Fs, A) as independent parameters. Here we 
define the origin of the s field to be the point where the messenger fields become massless, 
i.e., M mess = A(s). Once we integrate out the messenger fields, the gauge kinetic term in this 
case is given by 

'=-K?-(^ iog C) w '" w '" +h ' c ' <4) 

where g is the gauge coupling constant and N is the (effective) number of messenger fields. 
Our later discussion can apply when the low energy effective theory is of this type. For 



example, in the model of Ref. [35| the supergravity effects create a local minimum at (s) 



A 2 /Mpi, where Mp\ ~ 2.4 x 10 18 GeV is the reduced Planck scale. Ref. [36] discussed a model 
with an additional superpotential term, W 3 M mess ff, with which the effective value of (s) 
is M mess /A. 

The important point here is that the scalar field s couples to the MSSM fields with 
a suppression of Fs/(s} 2 , whereas the coupling to the gravitino is suppressed by F5/A 2 . 
Therefore, for (s) <C A, there is a possibility to avoid the dangerous gravitino overproduction 
as well as the catastrophic entropy release from the s decay [37, 27, [2^, |30|, [D], [l6|, 17, 31]. 

The three parameters (s), Fs, and A can be expressed in terms of physical quantities 
relevant for our discussion, such as the masses of Bino, s, and gravitino [m^, ms, ^3/2]- 
They are related to (s), Fs, and A as 

m~ - (5) 
mB ~ (4vr)2( s )' (5j 

m s = — , (6) 

m3/2 = tS^' (7) 

where g\ = ^5/3gy with gy being the coupling constant of the U(l)y gauge interaction. 
When expressed in terms of the running Bino mass at the electroweak scale, explicit A 
dependence of the low energy quantities disappears in many places. We can invert the above 
relations and write (s), Fs, and A in terms of the physical quantities: 

Fs = UxlO»G.V>( i ^), (8) 
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Here and in what follows, we use m 3 / 2 = 30 MeV, = 200 GeV and mj = 100 GeV as 
reference values, though the following discussion is generic and does not depend on those 
explicit values. 

Although A does not appear in the above relations, it cannot take an arbitrary value. 
In fact, there are lower and upper bounds on A to avoid instabilities at the SUSY breaking 
minimum. In order to avoid a tachyonic mass for the messenger fields, A should satisfy 
A 2 (s> 2 > XF S , i.e., 

A > A min = 1.7 x 10- 7 • iV- 2 ( m3/2 ) 1 ( m % ) 2 . (11) 

V30MeV/ V200 GeV/ V ' 

On the other hand, the interaction term in Eq. (0) induces a logarithmic potential at one-loop 



level H, HI, 

V(s) = ml\s - <*)| 2 + ^A 2 Fjlog (^) , (12) 

where we have assumed that the messenger fields transform as 5 and 5 under SU(5). The 
logarithmic potential gives an attractive force on the s field toward the SUSY vacuum at the 
origin. The stability at the SUSY breaking minimum requires 

5 V A 2 F| < \ml{ S )\ (13) 



167T 2 * 4 

namely 

A<A„^ 1 , Xl 0-3.^ ( _^ )( _^_)-. (14) 

In the following, we assume A m i n < A < A max . The above logarithmic potential also induces a 
mass splitting between the real and imaginary parts of the s field, 5m? /m| = 0(A 2 /A 2 nax ), as 
well as a shift of the minimum, 5(s)/(s) = C(A 2 /A max ). For simplicity, we assume A 2 <C A^^ 
and neglect those corrections in the following discussion. 



3 Scenario 



Let us first give an overview of the cosmological scenario in this model; (i) the s field develops 
a large expectation value during the inflation, (ii) its coherent oscillations after the inflation 
dominate the energy density of the Universe, and (iii) its decay produces radiation (including 
SUSY particles if kinematically allowed) and gravitinos. The assumption (i) is quite natural 
as far as ms is much smaller than the Hubble parameter H during the inflation, since the 
minimum of the potential during and after the inflation can be well separated from s = 
due to the deformation of the potential through gravitational (or general 1/Mpy suppressed) 
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Figure 1: A schematic figure of the evoiution of s field. 



interactions. The s field then starts to oscillate around the minimum when the Hubble 
parameter becomes comparable to ms and keeps oscillating until it decays. 

In the rest of this section, we discuss several conditions for the above scenario to work. 
An important point here is that there is a global SUSY minimum of the potential at s = 0, 
apart from the local SUSY breaking minimum, s = (s). Hereafter, we take the basis where 
(s) is real. (See Fig. |].) As discussed in Ref. [32], the s field does not fall into the SUSY 
vacuum unless its initial value Si n ; is too close to the real axis. We investigate in more detail 
the conditions for s to be trapped at the SUSY breaking minimum. When discussing the 
dynamics of the s field, we neglect the corrections from higher order terms in the Kahler 
potential in Eq. (|l]), which is small as far as \s\ < A. 



3.1 Avoiding the SUSY vacuum 

If the value of s approaches too close to s = during the oscillations, the scalar components 
of the messenger fields may become tachyonic, which makes the Universe quickly fall into the 
SUSY vacuum. This can be avoided if the initial value for Im[s] is so large that the trajectory 
of s stays away from s = 0, satisfying 



\s\ > 




(15) 



in the course of oscillations. 

Even if the above condition is met, the motion of s can be significantly affected by 
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the deformation of the potential near the origin due to the logarithmic potential Eq. (|12|). 
Moreover, it is known that a scalar field oscillating on a scalar potential of the logarithmic 
form experiences strong spatial instabilities and quickly deforms into spatially random and 
inhomogeneous state [pq] . If such instabilities becomes significant before the s field gets 
trapped in the SUSY breaking minimum, we expect that it falls into the SUSY vacuum. 
This can be avoided if the logarithmic correction remains subdominant along the trajectory 
passing near the origin, i.e., |s| <C (s), and the condition is given by 

16ir z m s 

More rigorous derivation of (16) can be found in Appendix [A|. 

The minimum value of |s| during the oscillations is approximately given by (cf. Fig. ||) 

,Im[s ini ] 

Mmin ^ (s)—. p, (17) 

I ^ini | 

where Si n ; is the value of s when it starts coherent oscillations at H ~ ms- We can rewrite 
the constraints (^) and ([!(]) respectively in terms of the ratio of the initial amplitudes, 
r = |Im[s in i]/Re[s in i]|, as 

A V 1/2 / ms \ / m 3/2 n -1/2 



v / TTr 
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* — I/2 (l^)(^v)(l^v)" 



In order to study the spatial instabilities, we separate the s field into a homogeneous part 
s and a perturbation <5s. We have numerically followed the evolution of s and 5s for a set 
of reference values of m 3 / 2 , wig, ms, and N, keeping only terms linear in 5s. The initial 
conditions are set as Re[si n j] = \/2Acos0 and Im[sj n i] = y^Asinfl, where is related to r as 
tan^ = r. We have chosen the initial value of \5s\ equal to 10~ 5 A, but the following result is 
not sensitive to this value. (Indeed, we have confirmed that the result remains almost intact 
for \5s\ = 10 _10 A.) In Fig. ^ we show a parameter region surrounded by a solid (green) where 
5s remains smaller than s until the s field settles down in the SUSY breaking minimum. The 
allowed regions are found to be 1(T 7 < A < 1(T 3 and 0.03 < 9 < vr/2. For different values 
of 7713/25 m B ano - m S' t ne allowed regions are modified correspondingly; in particular, the 
minimum value of 9 can be smaller. Note that 5s ~ s does not necessarily mean that the s 
field falls into the supersymmetric vacuum. What we would like to emphasize here is that 
there is a parameter space where our scenario is realized. We have also shown in the figure 
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Figure 2: The region surrounded by the solid (green) line represents a parameter space 
where the perturbation 5s remains smaller than the homogeneous part s until the s field gets 
stabilized at the SUSY breaking minimum. The conditions (18) and ( |l9| ) are also shown as the 



dot-dashed (red) and dashed (blue) lines. We have chosen m 3 / 2 = 30 MeV, = 200 GeV, 
m s = 100 GeV, and N = 1. 

the conditions (^) and fll9|), and the latter gives a slightly milder constraint on A than the 
solid (green) line. 

Another concern is whether the ratio of the energy densities of Im[s] and Re[s], r^ s = 
p Sl /p SR , is conserved or not. According to our numerical calculations, the final value of r c g is 
always larger or equal to r, and r e g- tends to become larger for smaller 9 and larger A. For most 
of the region surrounded by the solid (green) line, however, the ratio does not significantly 
evolve, and in particular, it remains almost constant in the course of evolution for 6 > 0.1 
and A < 10 -4 . Therefore we do not distinguish r e g- from r in the following discussion. In the 
above analysis we have not taken into consideration thermal effects, which will be discussed 



in Appendix B.2. 



3.2 s-dominated Universe 

We assume that the initial value of \s\ is so large that the coherent oscillations of s dominate 
over the energy density of the Universe before the time it decays. Such a domination happens 



8 




Figure 3: The condition of the s-dominated Universe. 

if 



2 



rp(R) - ■ frp rp X f \^[s ini ] -{s)\ \ 

T d < mm{T R , Tmg ){ ) (20) 



or 



T p < mm(rR , r „ s) (MM) 2 (21) 



where T^ 1 ' are the decay temperatures of sr/i (cf. next section), Tr is the reheating 
temperature after the inflaton, and T ms is the temperature at H = ms in the radiation 
dominated Universe, 

TTO , 7xl „» Gev( ^ ( _^)^ (22) 

where 5* counts the relativistic degrees of freedom in plasma. 

As we will see in the next section, T^p is always smaller than or equal to 2j . The 
history of the Universe depends on the values of the decay temperatures and r. There are 
four possibilities (cf Fig. ||): 

Case 1: [s/-domination] For r > 1, the energy density of sr is always smaller than 
that of si, irrespective of whether ( p0|) is satisfied or not. The si dominates the energy of 
the Universe, if ( ^T| ) is satisfied. 



Case 2: [s/-domination after SR-domination] If (|20|) is satisfied for \lT^ /T^ < 
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r < 1, it is sr that dominates the energy of the Universe first, and the si dominates the 
energy of the Universe after the decay of sr. 



Case 3: [^-domination] If (f20|) is satisfied for r < y /T^ R \< 1), sr dominates the 
energy of the Universe, while si does not. 



Case 4: [s/-domination] If (|20| ) is not satisfied while (|2l| ) is satisfied for r < 1, sj 
dominates the energy density of the Universe. 

Once the energy density of the Universe is dominated by either sr or s/, later discussions 
of the non-thermal gravitino production do not depend on the reheating temperature Tr. In 
Appendix [b], we will consider the thermal effects, and a consistent parameter region for the 
s domination is found to be 1(T 5 < A < 1(T 3 , 0.2 < 9< tt 72, and 10 5 GeV < T R < 10 6 GeV 
for the set of reference values: m 3/ / 2 = 30 MeV, m B = 200 GeV, ms = 100 GeV and N = 1. 
We would like to emphasize here that the consistent ranges for A and 9 depend on the choice 
of m 3 /2, fn B and mj. For instance, the lowest allowed value of 9 can be as small as 0(0.01). 

4 Decays of the s field 

The s field mainly decays into the MSSM particles through loop diagrams of the messenger 
fields. We first discuss the main decay mode and calculate the decay temperatures. 

4.1 Decays of s R 

The effective couplings between sr and the MSSM fields can be read off from the (s) 
dependencies of low-energy parameters. For scalar fields, the interaction terms are given 
by 



V2( 



(fh2 



£ -=^Mi - aRlfr (23) 

( f\ 

The effective mass parameter (m cS ) 2 is a part of the scalar mass that is proportional to 
1/| (s)\ 2 , i.e., (rrigg ) 2 = —dm 2 -/d\og \{s)\ 2 . If gauge mediation is the only contribution to the 

( f) 

scalar masses, rnj s is identical to their masses. In realistic models of gauge mediation, the \x 
parameter needs to be generated by some mechanism, and such contributions to the masses 
of the Higgs fields may be independent of (s). 

The couplings to the gauginos, A, are 

42 = -P- • - s -rAA, (24) 
int y/2(s) 2 H 5 v ' 
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where m\ is the gaugino mass (cf. Eq. (H)). There is a similar coupling between sr to 
Higgsinos: 

£ (m ggsi „o) = --^. SR (W d . P L h u ) + h.c. (25) 

The coefficient /j, e g is again a part of \x that is proportional to l/(s). 

There are couplings to the quarks and leptons through a mixing between sr and Higgs 
bosons (h° and H°). The mixing is induced through the interaction term in Eq. (p3|) for 
/ = H u ,Hd with one of the Higgs fields replaced by its vacuum expectation value. 

The couplings to the gauge bosons are 

{A) _ 1 2g\N 1 1 

Ant " V^J^fV) ' a SrF ^) f ^ (26) 

where the index A represents the gauge group (SU(3), SU(2), and U(l)). 

The main decay mode of sr depends on the mass spectrum of SUSY particles. We discuss 
the case of the sweet spot SUSY model |$^] in detail as an example of realistic models. In this 
model, the ^-parameter and the Higgs soft masses are generated at the GUT scale through 
direct couplings between the Higgs fields and the SUSY breaking sector. Those contributions 
do not depend on (s). Additional large contributions to the H u soft mass are generated 
through gauge mediation and the renormalization group (RG) running. In particular, there 
is a significant RG effect due to the large Yukawa coupling of the top quark and the large 
scalar top masses. These contributions are proportional to l/|(s)| 2 , and thus enhance the 
effective coupling to sr. The effective mass parameter ("i^ff"^) 2 is estimated to be 

- {m^ f = K) 2 , (27) 

with 

k ~ 3 - 4. (28) 
The parameter k depends logarithmically on the messenger scale. For the down type Higgs, 

\(rnS d) ) 2 \<<\(m^n (29) 
due to relatively small RG effects. The effective coupling to the Higgsinos is also suppressed, 

l^effl < N> (3°) 

in this model. 
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The enhancement in Eq. (pg|) is very important since there is no such factor in the sj 
decay. The decay modes sr — > hh, ZZ, and WW (where the gauge bosons are longitudinally 
polarized) and also the fermion modes such as sr — > bb through the SR-h mixing are 
enhanced. This makes the sr decay much faster than that of sj in the parameter region 
of our interest. 



The partial decay width of the hh + WW + ZZ mode is given by 

i ( V2(m3H; 

8Trms \ (s) 



1 (V2(m ( ^ ) ) 2 S m 2 p\ . . 

s R -+hh + Ts R -^ZZ + l sr -+WW — o 7~\ • 



We ignored the mass differences among h, Z and W, and also C(m|/m|) and 0(w? h jrn?^) 
terms for simplicity. The angle j3 is defined by tan/3 = (H u ) j {H^) , and tua is the mass of the 
pseudo-scalar Higgs boson. If sr mainly decays into hh + WW + ZZ, the decay temperature 
T^' is given by 

T iR) ~ 13GcV Ar-i^g^' 1 / 4 / m e? u) rn s yi/y m 3/2 



d ' V15/ ^800GeV i V500GeV/ V30MeV 

(32) 



,200 GeV/ 

where we have defined the decay temperature as[j] 

7f } - (^)^ /4 7rVV%. (33) 

In the numerical calculations below, we have included the phase space factors and contribu- 
tions from other decay modes such as ti. 

The width of the bb mode through the mixing to the lightest Higgs boson is given by 

3m s / Vgmg^g sin 2 (3 
8n { (s) ml 

Correction terms of C(m 2 /m|) are ignored. When the bb mode is dominant, the decay 
temperature is given by 

rf) * i.6 Gev • iv- 1 ( *r 1/4 ( V ( mh y 2 ( ms v /2 

d U5/ UoOGeVj V115GeW VIOOGeV 




m 3/2 x-i 



30MeV/ V200GeV 



2 

1 £ 

2 



sin 2 /3. (35) 



^Note that this is a temperature of the Universe when the age of the Universe is comparable to the lifetime 
of s R , provided that the Universe is radiation dominated or radiation produced by the s R decay dominates over 
the Universe. Although does not represent a temperature if the decays happen during the si dominated 
era, we call this quantity the decay temperature in later discussion even in that case. 
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3 N = 1 tan/3 = 35 \m eS \ = m A = 4m^ 3 N = 1 tan/3 = 35 l^eff I = m A = &nig 




m s [GeV] m s [GeV] 

Figure 4: The decay temperatures of sr and s/. We fixed the gravitino mass to be 30 MeV. For 
other values of m^j 2 , the decay temperatures are obtained by multiplying (m 3 / 2 /30 MeV) -1 . 

We show in the left panel of Fig. Q the decay temperature as a function of ms 

for nig = 100, 200, 300 GeV. In the figure, we fixed the gravitino mass to be 30 MeV. 

The decay temperatures for other values of m 3 / 2 can be obtained by multiplying a factor of 

(7713/2/30 MeV) -1 . We have fixed other parameters to be N = 1, tan/3 = 35, nth = 115 GeV 
(H ) 

and |»rig£f I = 777^4 = Am B , where 777,4 is the pseudo-scalar Higgs boson mass. In the calcula- 
tion, we have included the decay modes sr — ► hh, ZZ, WW, 77, BB, WW, gg, ti, bb, tt, ff, gg. 
The main decay mode of sr is sr — ► hh, WW, ZZ for 2mw < ms < 1 TeV. The gaugino 
modes become important for 7775 > 1 TeV. For 7775 < 2mw, the sr — > bb decay through the 
SR-h mixing is the main decay process. We can see a sharp peak at 7775 = mh due to the 
enhancement of the mixing. 

4.2 Decays of sj 

In Ref. |32|, the decay property of si is assumed to be the same as the one of sr. We show in 
this subsection that the decay of s/ happens much later in particular when sr — > hh is open. 
The difference of the decay temperatures will be important in calculating the non-thermal 
gravitino abundance. 

The imaginary part, s/, can only couple to CP-odd combinations. The enhanced coupling 

( H } 

to Higgs bosons through a large value of m eS u is therefore absent. There is a coupling to 
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the Higgs bosons through an (s) dependence of the Bfi-term: 

i m \ sin 23 

C 3 BfiH u H d + h.c, -► (s) SlHuIId + hx ' (36) 

Here we have assumed Bfi oc (s)" 1 and used a tree-level relation from electroweak symmetry 
breaking, Bfi = — (m^sin2/3)/2. The coupling constant is suppressed by a sin 2(3 factor 
which is generically small in gauge-mediation models. 

The decay into two gauginos is therefore important if it is kinematically open. The 
interaction Lagrangian is given by 

Ant 3 ^J^' 2 Sl ^ ilhX - ( 37 ) 

Since this is the same strength as the sr\\ coupling, the lifetime of si is always longer than 
sr due to the suppression of the Higgs modes. The partial decay width of the Bino mode is 
given by 

r _m,(m B \'( 4,„|y' 2 

r ».-^- 32^77 ) y-^i) ■ (38) 

If this is the dominant decay channel, the decay temperature is given by 

l d _ uuMeV iv ^ 500 GeV J UoMeVv 1 V200 GeY J 

( 4m|\ 1/4 



where the decay temperature is defined in the same way as (p3[). The Binos subsequently 
decay into staus if > nif. In a large ms region, the Wino and the gluino modes become 
more important. However, as we will see later, such a region is not allowed because of the 
overproduction of gravitinos. 

If the Bino mode is closed, the main decay mode is into bb through the sj-A° mixing from 



Eq. (36). The partial width is calculated to be 



_ 3ms ( m\ sin 2 (3 m b . 
s '^ bb 16vr V <*> ' m\-m\ ' ' 



When the bb mode is dominant, the decay temperature is 

if - i6 Mev • at- 1 r*r 1/4 ( ms Y /2 ( m ^ rV m * 

d \lh) VlOOGeV/ V30MeV/ V200 GeV 



-i 



sin 2 3. (41) 
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The value in front becomes 18 MeV once we include the tt mode. 

The decay temperature as a function of ms is shown in the right panel of Fig. We have 
used a set of parameters which are indicated in the figure. Again, for other values of m.3/2, 
T\ oc m~j T As we can see, ^ is significantly lower than for mg < 1 TeVj^] The decay 
modes si —* BB, WW, gg, tt, bb, tt, gg, 77, hA, HZ are included in the calculation. We have 
ignored the hZ mode because it is much smaller than the bb mode. The gauge invariance 
requires the hZ mode to vanish in the decoupling limit, m^/m^ — » 0. 



5 Non-thermal gravitino production 

The s field can decay into two gravitinos with a suppressed branching fraction. We calculate 
here the branching ratio and estimate the gravitino energy density. We will see that the non- 
thermal component can explain the DM abundance when mj ~ 0(100) GeV independent of 
the gravitino mass. 



5.1 Abundance 

The non-thermal gravitino abundance can be calculated from the decay temperatures and 
the branching ratios of the srj — ► decays. The partial decay width of srj into two 

gravitinos, srj ^3/2^3/2, is given by ||, ||, [l7|, |2| 



r 3/ 2 = — V- ( 42 ) 

3/2 967TJW* \m 3/2 J ( ' 



This formula is obtained from the interaction term in the second term of Eq. (JTJ) by identifying 
the fermion component of S with the longitudinal mode of the gravitino. By using this partial 
decay width, we can calculate the branching fraction. 

There are two interesting branches where the main decay modes are different. For the 
sr decay, the main decay mode is (A) sr — ► bb for 2m;, < mg < 2niy/, and (B) sr — > 
hh, WW, ZZ for 2mw < m s i$ 1 TeV. The branching ratios of the two-gravitino mode in 
those cases are respectively given by 



b^ ~ 46xio- 9 -iv 2 f |m ^" }| 1 r mg vr m * r 2 f mfe v 

^3/2 - ^°xiu 800 GeV V100 GeV/ V 200 GeV/ Vll5GeV/ 




'"Two decay temperatures become similar when the gg mode becomes dominant, i.e., ms < 2m T . 
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< 2 ) * 2.2 x 10- 7 • iV 2 ( " f mg ) 6 f m * ) " 2 sin" 4 /? • • • (B). 
3 / 2 ^800 GeV^ V500 GeV/ V200 GeW M V ; 

It is interesting to notice that the branching ratio is independent of 1713/2- For the sj decay, 
the main decay mode is (C) si — ► 66 (ms < 2m^) or (D) si — ► -B-B (ms > 2m The 
branching ratios in two cases are 

Bg> - 4.2 x 1(T 5 • iV 2 f mg V ( ^ ) ~ 2 ( 1 - 4V sin" 4 •••(C), 
3 / 2 VlOOGeW V200 GeW V ™a/ 



2 \ -1/2 



We here define quantities fij^ and which represent the density parameters of the 



gravitino when we ignore the presence of sj and s^, respectively: 

d 

m s 



0$ ee -m 3/2 ^B<;»x2/w S ) 0l (43) 



n s/2 = 4 m 3/2^Bi; ) 2 x2/(p c / S )„, (44) 



where (p c /s)o — 1-8 x 10 9 GeV is the critical density divided by the entropy density at 
present. The abundances and fi|w 2 are related as 

= ^- (45) 

o (-0 ' 

12 3/2 

where we used the fact that B^/B^l = (rfp /T^) 2 . In the actual situation, of course, one 
cannot totally neglect sr or s/, and one has to take into account both of the contributions 
and also the dilution effects. The gravitino abundance in a general case can be expressed in 
terms of Sl{ R l and Q^L as 



"3/2 aiLKX "3/2 



i4 3/2 



r 



(46) 



The former and latter regions of r respectively correspond to the cases where sj does not and 
does dominate the energy density of the Universe before the sj decays. For r > T^p /T^ R \ 
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Figure 5: The non-thermal gravitino abundance with respect to ms for several values of the 
Bino mass = 100,200, and 300 GeV (Left). We have varied r from 0.01 to 0(1) for each 
value of rrijj. The other parameters are shown above the panel. For r > 1, the three regions 
overlap one another around ms ~ 100 GeV and > 100 GeV. The right panel shows the 
contours of = ^-2 for various values of r. For r > O(0. 1), the abundance does not 
depend on r since the sj oscillation dominates over the Universe. The gravitino abundance 
becomes larger as ms increases (see Eq. (fT7|)). 

most of the gravitinos are produced by the si decay. For r 2 > Tj 1 ' /Tj , both radiation and 
gravitinos arise from the sj decay, and thus the gravitino abundance becomes insensitive to 
r. Note that (^3/2) corres P on ds to the gravitino density parameter in the limit of r — > 
(r — > 00). 

Since T^' 1 ' oc l/(s) oc m~j 2 and the branching fractions B^j^ are independent of m 3 / 2 , 
both an d ^3/2 are a ^ so independent of the gravitino mass. Therefore, interestingly, the 
total gravitino energy density ti^J 2 does not depend on the gravitino mass. 

We show in the left panel of Fig. || the gravitino abundance for fixed values of 
for r > 0.01. The non-thermal gravitino can account for the observed DM abundance for 
ms = O(100) GeV. In the right panel, contours of ^^J 2 = 0.2 for various values of r are shown 
on the ms-m B plane. For r > O(0. 1), the contour gets independent of r for m^ > 100 GeV. 
For r > 0.05, the correct abundance is obtained for ms ~ 100 GeV where si — ► bb is the 
dominant decay mode. The gravitino abundance for that case is given by 

NT / mQ \7/2 / m ~ \ -l 



The abundance does not depend on the detailed model parameters such as k, mA or tan/? as 
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long as rriA 3> m$ and tan (5 > 3. In general, this result applies to models where the 5^-term 
is proportional to and the decay width in Eq. is larger than or comparable to that 



in Eq. (gD 



Although the abundance is independent of m 3 / 2 , the gravitino mass cannot be arbitrarily 
large. In order to avoid overproduction of 4 He, the decay temperature of si is required to 
be higher than ~ 10 MeV |39|], which gives an upper bound on the gravitino mass. For 
ms = 100 GeV and > 100 GeV with r > 0.05 motivated by DM abundance, we obtain 
from Eq. (|l|): 

™ 3 /2 ~ °( 10 °) MeV - ( 48 ) 

If one allows a small value of r, there is a region where sj — > BB is open while Ufjf = 0.2 
is satisfied. The upper bound on the gravitino mass in this case is relaxed to 0(1) GeV. Such 
a region is subject to the BBN constraint from the decay of the non-thermally produced 
NLSP. We will discuss the constraint later in the Appendix [C]. 

5.2 Thermal component 

Here we comment on the amount of thermally produced gravitinos. After accounting for 
the dilution effect by the entropy production of the s decays, the density parameter of the 
thermally produced gravitinos is estimated to be (3^j 

-2 



^ 3 h /2 ~ a004 ( 30 Mev) (200 Gev) ( 18 MeV ) {2.5 x IIP & GeV 1 ' (49) 

where |si n ;| is the size of the initial amplitude. If the entropy of the Universe is generated by 
the decay of sr (s/), one should substitute (T^) for T^. This expression is independent 
of the reheating temperature after inflation even though most of the gravitinos are produced 
at the end of the reheating process]^] 

Let us assume that s/ dominates the energy density of the Universe. In the case where 
the si — > bb mode is the dominant decay process, we obtain 

I \ 2 

L ^ " UW± ' iV VlWGeVV v3¥MeVV v2WGeVV IT 



where A is the cut-off scale in Eq. (|T]). The |si n i|/A factor cannot exceed 0(1) for the discussion 
to be within the framework of the effective theory. In order for the thermal component not 



§ There is a logarithmic dependence on the reheating temperature through the running coupling. 
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to exceed the observed DM density of the Universe, we obtain a lower limit of the gravitino 
mass: 



™3/2 > O(10) MeV, (51) 
for m s ~ 100 GeV and ~ 200 GeV. 



On the other hand, if sj — * BB mode is open, the gravitino abundance becomes 
'3/2 - U -- L - JV V5Q0 GeV J vlOO MeY J ^200 GeV, 





A ' 1 ' -» ■ ' < 52 > 
The gravitino mass must be slightly heavier than the previous case, 

m 3/2 > O(100) MeV, (53) 
for m s ~ 500 GeV and ~ 200 GeV. 

5.3 Free-streaming scale 

We have seen that the gravitinos non-thermally produced by the decay of s successfully 
account for the observed abundance of DM. Since the gravitinos are relativistic at the 
production, we need to check if the free-streaming scale is consistent with the observational 
bound, \ps < 0(100) kpc, from the Lyman a forest data p(|. 

Let us first derive an expression for the free-steaming length, A^s, which is a distance 
that particles (gravitinos in our case) can travel until they become non-relativistic. In the 
following the gravitinos are assumed to be produced by the decay of s (either sr or sj), and 
the decay temperature denotes either T^' or Tfjp . Let us denote by p the momentum of 
the gravitino. Due to the cosmic expansion, the momentum red-shifts as p oc a , where a 
denotes the scale factor. Noting that the velocity of a particle is given by the ratio of the 
momentum to the energy, the free-streaming length is expressed as 

where pd is the initial momentum at the production, ad and a eq are the scale factors at 
the decay of s and at the matter-radiation equality, respectively, and the scale factor is 
normalized to be unity at present. We take the matter-radiation equality as the end point of 
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integration, assuming that the gravitino has become already non-relativistic at the equality. 
The assumption is satisfied for the whole parameter space of interest. 

Assuming that the Universe was radiation dominated since the gravitino production until 
the equality, we can perform the integration and obtain 



1 + 



H, 



z eq ( m s a d 
V2m 3/2 a, 



eq 



~ 60kpc(|^ 



sinh 



m 



3/2 



30MeV 



2m 3/2 CLeg 

m s a d 

m s 
100 GeV 



1 -0.1 In 



™3/2 

30MeV 



-i 



18MeV 

ms 



100 GeV 



18MeVy 



(55) 



where z eg and H eq are the red-shift and the Hubble parameter at the equality, respectively, 
and we have adopted an approximation, ms 3> wi 3/ / 2 . In the last equality, we have used 



eq 



3176 and H eq ~ 31 Mpc" 1 



If the decay of s produces not only the gravitinos but also the (almost) entire entropy of 
the Universe, we can express the free-streaming length in terms of the branching fraction of 
the gravitino production. The free-streaming length is then give by 



60kpc 



x < 



1-0.1 In — 




(56) 



Since -B3/2 is independent of 7713/2, the free-streaming scale is also independent of 7773/2 once we 
fix the gravitino abundance. Note that the expression Eq. (|56|) is valid only for r < /T^ 



or r > yT^/T^ R \ while Eq. (55) holds for any values of r as long as the gravitino comes 
mainly from either sr or s/. In the parameter region of our interest, where si — > bb is the 
main decay mode, the free-streaming scale is of 0(100) kpc, which is on the border of the 
Lyman-a bound, Xfs 0(100) kpc. It is an interesting possibility that we may be able to 
see the suppression of the structure formation below the corresponding scale. 



5.4 Isocurvature perturbations 

During inflation, s is assumed to be at |sj n i| ~ A, far deviated from the origin. If the s field 
has an approximate U(l) symmetry, the phase component, 6 = arg[s], remains light and 
therefore acquires quantum fluctuations 56 = il//(27r|si n i|), where Hi represents the Hubble 
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parameter during inflation. The phase 9 is related to r (the ratio of the initial values of sr 
and si) as tan# = r. Therefore, 59 amounts to the fluctuation 6r, which generically leads to 
the isocurvature fluctuations in the gravitino DM. This is because sr and s/ have different 
decay temperatures and different branching ratios into the gravitinos. 



Recall that the gravitino abundance becomes insensitive to r for r > y /T^ R \ This is 
because both the radiation and the gravitino are produced mainly from the decay of si. Thus, 
the cold DM (CDM) isocurvature perturbation is also suppressed in this case pi] . Intuitively 
speaking, for a large enough value of r, we can simply neglect the sr; the fluctuations in 
radiation and the gravitino DM are then adiabatic, since both are generated from a single 
source, sj. 

To see this more explicitly, let us estimate the CDM isocurvature perturbation S c ~, in the 



case of r > ylT- f'/T\ } . From Eqs. © and (46), we have 



where we have used 5r = 59/ cos 2 9 in the last equality. As we mentioned above, we can 
see that S n is suppressed for r ^> \l /T^ R \ The current observation bound on the 
isocurvature perturbation reads (.S^l ^2x 1CP 5 at 95%C.L. Q. Thus, Hj is bounded above 
not to exceed the current constraint on the isocurvature perturbation, 

iU .r<.*r( r3 \ ( T d R) \ ( m 3/2 W m s \ -l 



El < 2 x 10 " gov [j^)\^m) (sow ) lioo^vJ IT J ■ < 58 > 

There are plentiful inflation models which satisfy the bound. 

6 Sweet spot 

We have seen that the non-thermal production of the gravitino can explain DM of the Universe 
in a class of gauge mediation models. For r > 0.05, the correct abundance is obtained when 
the si — > bb mode is the dominant decay mode. The abundance in that case is given by 
Eq. (|47| ) which is independent of 7713/2- The constraints from the thermal production of the 



gravitino (Eq. (51)) and the decay temperature (Eq. (|4q)) restrict the mass range of the 
gravitino to be 

10 MeV < m 3/2 < 100 MeV. (59) 
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Interestingly, the above mass range and mj ~ 100 GeV overlaps with the prediction of the 
gravitational stabilization mechanism in Ref. p5[. This model relates (s) and A by 



W = (60) 
x 1 6Mpi v ' 

which is translated into a relation among ms, w-3/2 an d as 

The reference value we took approximately satisfies the relation. It is also interesting to note 
that the above supergravity effect always exists. Therefore, there is no big room left for other 
mechanisms to give messenger masses in the scenario of the s-dominated Universe. 

Ref. |33| proposed a solution to the [i problem by using the above gravitational stabi- 
lization mechanism. The //-term is generated from the direct interaction terms between the 
SUSY breaking sector and the Higgs fields, K 3 S' H u Hd/ A. This framework predicts 

^ ~ X ~ ms ' ^ 

which is perfectly consistent with fx ~ 0(100) GeV required from electroweak symmetry 
breaking and ms ~ 0(100) GeV from gravitino DM. 
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A Spatial instabilities 

Let us estimate a condition for an instability to grow. The scalar potential of s = (sr + 
isi)j\pl is given by 

V(s R , sj) = \ml ((s R - V~2{s)f + sj) + ^A 2 Fj log (%^) • (63) 
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Differentiating the scalar potential with respect to sr and sj, we obtain 
u -z _ 1 7)3^ ds R d Sl 

M = I d 2$ Q2y 

dsj{dsj ds'j 

m 2 | 57V \2 p2 ( 2 4s fl \ 5iV x 2 F 2 4s fl 3j 



5W X 2 F 2 4sflSj 2 , 5iV X 2 F 2 / 2 4s ? 



(64) 



Neglecting the cosmic expansion, the instability grows if det[.M 2 ] < 0. Thus, the condition 
for the instabilities not to grow is det[7W 2 ] > 0, namely 



16^ 2 m| ft \\ mSiX J 4 ' 
or equivalently, 

A < 2A max sin0. (66) 



B Remarks on initial conditions 



In this Appendix, we discuss conditions for srj to dominate the energy density of the 
Universe, taking account of finite temperature effects. For concreteness, we set m.3/2 — 
30 MeV, nig ~ 200 GeV, ms ^ 100 GeV and N = 1, and also r > 0(1) as reference values, 
which lead to a successful gravitino DM scenario from the si decay, as discussed in the text. 
Therefore we are concerned with a condition for si to dominate the energy density of the 
Universe. 



B.l s/-domination 

The ^-domination condition Eq. ([H]) leads to 

|ImM > 7xl „ 12GeV .(=&i)- i/2 (_^_y /2 . (67) 

If we take a natural expectation, |si m | ~ A, the si domination can be realized for Tr > 
10 5 GeV. 



B.2 Finite temperature effects 

In this subsection, we assume that the s field starts its oscillations when the Universe 
is dominated by the oscillating inflaton. Even before the reheating, however, there is a 
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background dilute plasma with a temperature T ~ (T^Mpi-ff) 1 ' 4 . The potential of s field 
therefore receives thermal corrections, which are not taken into consideration so far. Here we 
briefly discuss the finite temperature effects on the evolution of the s field and the messenger 
fields. 

There are two thermal effects on the s field: thermal mass and thermal logarithmic terms, 
which arise depending on whether the messenger fields are in thermal bath or not. If the 
effective masses of the messenger fields are smaller than the temperature of thermal plasma, 
i.e., m e ff = A|s| < T, the messenger fields will be in thermal equilibrium. The s field then 
receives a thermal mass: 

V± l) = ^iVA 2 ^!^) 2 for A|s| < T, (68) 

where we have assumed that the messenger fields transform as 5 and 5 under SU(5). 

On the other hand, when the messenger fields are so heavy that they are decoupled from 
thermal bath, there is a thermal effect arising from the two-loop contribution to the free 
energy, 5V oc g(T) 2 T 4 . Here we consider only the SU(3)c gauge group, which gives the 



dominant contribution to the free energy [42] 



5V = yflg(T)T*. (69) 
For A | s | > T, the running gauge coupling g% (T) is modified as 

9s(T)\x\ s \>T = 53(T)| s=(s) +iV^0l!ln(^), (70) 

where Mjj is some ultraviolet scale where g% is fixed. This leads to a thermal correction to 
the scalar potential 

V? = ^«3(T) 2 T 4 ln(Mj) forA| S |>T, (71) 

which may become important where the thermal mass term is negligible. 

Next let us consider the thermal effect on the messenger fields. The messenger fields 
acquire thermal masses through the gauge interactions with the SSM particles in thermal 
plasma. The thermal masses tend to prevent the messengers from falling into a SUSY 
minimum. In principle this effect could enlarge the allowed region for r: even with a small 
value of r, the messengers may be stabilized at their origin and the s field may settle down 
at the SUSY breaking minimum in the end. However, if this is the case, our scenario would 
be modified in two ways. First, the messenger fields are in thermal equilibrium when they 
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are stabilized by their thermal masses. If the messenger number is conserved, the lightest 
messenger may exceed the DM abundance. Although this issue can be avoided by introducing 
the breaking of the messenger number, it would make the analysis model-dependent. Second, 
the gravitino abundance ( f49[) is modified because the gravitinos are also generated from the 



scattering processes including messengers [43|. Thus, in order to keep the success of our 
scenario in the text, we assume that the messenger fields are so heavy that they are always 
decoupled from thermal plasma. 

We have shown in Fig. || a parameter space in which (i) the perturbation 5s remains 
small compared to s until the s field is stabilized at the SUSY breaking minimum and (ii) 
the messengers remain decoupled from thermal bath during the course of evolution, for Tr = 
10 5 GeV and 10 6 GeV with the thermal effects (^) and (|7l| ) taken into account. Here we have 
set 7713/2 = 30 MeV, = 200 GeV, ms = 100 GeV and N = 1. Notice that, compared to the 
zero-temperature result shown in Fig. g, smaller values of A are excluded since the messengers 
would be thermalized. The allowed region disappears for Tr > 10 7 GeV. Therefore, our 
scenario works for 10~ 5 < A < 10~ 3 , 0.2 < 9 < vr/2 and 10 5 GeV <T R < 10 6 GeV, if thermal 
effects are taken into account. Note that the consistent ranges for A and 9 depend on the 
choice of ^3/21 m B an d m s- F° r instance, the lowest allowed value of 9 can be as small as 
0(0.01) for e.g. m s = 400 GeV and = 100 GeV. 



C BBN constraints 



If the decay of sr or si into the superparticles are kinematically allowed, they are copiously 
produced, which will decay into the NLSP promptly. Depending on the lifetime of the NLSP, 
their abundance is subject to the BBN constraint. We assume in the following that the NLSP 
is the stau since the constraint is much weaker than the Bino NLSP case. 

Before proceeding, let us mention when the BBN constraint could become important. As 
one can see from the right panel in Fig. ||, there are parameter regions where the decay into 
the superparticles is significant for rC 1 while the gravitino abundance is fixed. On the 
other hand, for the reference values of ITI3/2 = 30 MeV, = 200 GeV, ms = 100 GeV and 
iV = 1, we have found that r must be larger than 0.2 (0.05), if thermal effects are (not) taken 
into consideration for our scenario to work (see Figs. [2] and ||). However, for a different choice 
of those parameters, the smallest value of r can be O(0. 01). Thus, for a certain fraction of 
the parameter space of our concern, the BBN constraint may be important. 

The main source of stau is the sj — > BB decay followed by the decays of Binos into staus, 



25 




Figure 6: The parameter region where the s field is trapped by the SUSY breaking minimum 



and the messenger fields are decoupled from the thermal bath, when the thermal effects (68) 
and ( |7l| ) are taken into account. We have set m 3 / 2 = 30 MeV, = 200 GeV, mg = 100 GeV 
and N = 1. 

if the Bino mode is open. In this case, the decay of sr is not important since the si decays 
much later. If the Bino mode is closed, the main source is sr — ► ff. 

By using the decay temperatures calculated before, the non-thermal stau abundance can 
be estimated through (82) in Appendix O, where a general formula of the non-thermal relic 



abundance is derived. When a large number of staus are produced by the si decay, the fast 
pair annihilation processes make the final abundance approach a value determined by the 
decay temperature and the annihilation cross section, which is not sensitive to the initial 
abundance. 

We show in the left panel of Fig. ^ the abundance of the non-thermally produced stau 
as a function of mj with the same set of parameters as Fig. ||. The right figure is the case 
with m 3 /2 = 100 MeV. We have used the annihilation cross section of the staus in Ref. 4C]. 
(Recently it has been shown that the cross section can be larger if there is a significant left- 



right mixing in the stau sector [47, |48|].) The parameter r, the ratio of the amplitudes, is 
taken to be r = 0.01. 

The abundance Yf does not depend on r if the si — > BB decay is kinematically allowed 
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m s [GeV] m s [GeV] 



Figure 7: The abundance of the non-thermally produced staus. We took m? = 100 GeV. The 
gravitino mass is m 3/ / 2 = 30 MeV (left) and 100 MeV (right). The abundance is inversely 
proportional to the decay temperature. 

(ms > 2m,g). In the case where sr — > f f is the main production process (2mf < ms < 2m g), 
we should take into account the entropy production from the si decay which happens at a 
lower temperature. Therefore, for a larger value of r, the stau abundance in the ms < 2m g 
region is more suppressed by a larger dilution effect. For different values of the gravitino 
mass, the stau abundance approximately scales as Yf oc T7" 1 oc m 3 / 2 . 

For the parameter set we took, mf = 100 GeV and m 3/ / 2 = 30 MeV, the staus decay 
rather early in the BBN era (jf = 50 sec), and thus there is no significant constraint on 
Yf from BBN. The bound that the gravitinos from the stau decays should not exceed the 
observed matter energy density gives If < 1 x 10~ 8 , which is satisfied for any value of m$- 
For m 3 /2 = 100 MeV and mf = 100 GeV, the stau lifetime is 600 sec, with which we obtain 



a BBN constraint from the D abundance, Yf < 1 x 10~ 10 45]. In this case, a part of the 
parameter region is excluded as shown in the right panel of Fig. [7[ For a further large 
value of m 3 /2, the constraint from the 6 Li abundance becomes important. For example, for 
m 3/2 = 300 MeV, the constraint is Yf < 1 x 10~ 13 [fHf ], and the consistent parameter region 
disappears for mf = 100 GeV. 
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D Non-thermal relic abundance 



In this appendix we calculate the non-thermal relic abundance of a particle X, assuming 
the following cosmological scenario, (i) The energy density of the Universe is dominated by 
a non-relativistic matter (j) (e.g. a coherently oscillating scalar field), (ii) The eft field then 
decays into radiation and X, with a decay rate (iii) The subsequent pair annihilations 
of the X particles reduce its number until it freezes out. The relevant Boltzmann equations 



are given by (cf. |Tj, p0| ) 

Ptf, = -SHp^ - r^p^ , (72) 

Prad = -4-ffprad + ^<j>P(j> , (73) 

h x = -3Hn x - (crv)(n 2 x - n 2 x Jj + T^—b, (74) 

H 2 = 2 Ptotal; Ptotal = P4> + Prad , (75) 

6M p\ 

where p^ and p ra d are the energy density of the <j> and radiation, respectively, and we assume 
that the energy density of the X particle is negligible compared to them, px <C ptotal = 
Prad + P<j>- H is the Hubble parameter, nx is the number density of X, and b is the averaged 
number of X particles produced per <j>. Here and in what follows, we assume that the 
equilibrium number density is negligible, nx,eq *C nx, which is a good approximation as long 
as mx 3> Td with a moderate value of where is the decay temperature, defined by 
Td = (7r 2 g*/90) -1 / 4 ^jMpiT^. In terms of the following variables, 



the equations (|72|)-(|7q) become 



l- f i)^ = -\ U e* + l - U{ l-U), (77) 
1 f A dNx 3 f N N * +Af W 2 (7R\ 



A J dx 8 J(p 2e x / 2 

This can be solved numerically with initial conditions f<j>{— oo) = 1 and N x {—oo) = 0, and 
the final answer Nx(oo) depends only on the dimensionless parameter A, which is given by 

3M 2 r^(<™) 



A = 



2mj 



g*\i ( T d \ 2 ( (av) \/500 GeV 



77 * I0 ° ( SJ ' UowJ iiF^J l^) (A) ■ (79) 



'Note that this is different from the case of Q-ball decay , in which the X production suddenly terminates 
at T = Ta- The final X abundance obtained here is about 5 times larger than the case of Q-ball decay. 
"Strictly speaking, b must satisfy b ^> (m^m 3 ^ 2 /T^ 2 ) exp(—mx /Td). 



28 



The dependence of Nx(oo) on A is actually very weak, and it is empirically found that 

AT x (oo) -4.5 (l + 0.043 log (7^)) ■ (80) 

This approximation reproduces the numerical result within a few %, for a wide range of 
A = 10 3 — 10 10 . Thus, the final X abundance (for x — > 00, — > 0) is given by, assuming 
g* = const., 

n x A7 H 3 / 2 __ f 45 \ 1/2 iVx 



^rf^), Uvr^J M p iT(i(av) (81) 



~ 4.4xl0- n (-l ( ' ) ( - ' ) ( 1 + 0.043 log (-^1 1 (82) 



/10\ 1/2 fWO MeV\ flQ- 7 GeV~ 2 \ ( ( A \ 
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